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Cross sections for e + e~ — » hadrons, ir + ir~ J/tp, and fi + fj,~ have been measured in the vicinity of the ip(2S) 
resonance using the BESII detector operated at the BEPC. The ip(2S) total width; partial widths to hadrons, 
ir + n~ J/ip, leptons; and corresponding branching fractions have been determined to be V t = 264 ± 27 keV; 
F h = 258 ± 26 keV, T„+ v - J/4 = 85.4 ± 8.7 keV, and T t = 2.44 ± 0.21 keV; and B h = (97.79 ± 0.15)%, 
B^+n-j/i, = (32.3 ± 1.4)%, and Bi = (0.93 ± 0.08)%, respectively. 
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1. Introduction 

Since the discovery of the ip(2S) in 1974 @], 
measurements of its total width (T t ), and par- 
tial decay widths into hadrons (I\), 7r + 7r~J/^ 
(T 7T + 7T - ;u), and l + l~ (Ti), and the correspond- 
ing branching fractions, Bh, B^+^-j/^, and Bi, 
have been carried out [2—10]. The results of 
these experiments differ on both decay widths and 
branching fractions. The parameters are of par- 
ticular interest because, for instance, '0(25') — > 
is used in reconstructing B mesons for 
CP violation measurements pfl , and ip(2S) — > 
ir + ir~J/ip is often used to determine the total 
number of ip{2S) events in tj){2S) branching frac- 
tion measurements due to its large branching frac- 
tion and straightforward detection. Therefore, it 
is important to measure these decay widths and 
branching fractions with better accuracy. 

Twenty four center-of-mass energy points were 
scanned in the vicinity of the ift(2S) peak rang- 
ing from 3.67 GeV to 3.71 GeV. The data were 
collected with the BESII (BEijing Spectrometer) 
detector at the BEPC (Beijing Electron Positron 
Collider) storage ring. The BESII detector is 
described in detail in Ref. fL2| , In addition, 
separated-beam data were taken at the first and 
the last points for background studies. The total 
integrated luminosity was 1149 nb _1 . 

2. Event selection 

The following four reactions are studied: 

e + e~ — » e + e~ , 

e + e~ — > fi + \i~ , 

e + e~ — > ir + ir~J/i[i , 

e + e~ — > hadrons 

For the selection of lepton-pair final states, two 
charged tracks with total charge zero are required. 
For n + /J.~ events, the acollinearity must be less 
than 10 degrees. In addition, in order to suppress 
cosmic ray background, the time-of-flight mea- 
surements of the two muon-candidates must sat- 
isfy "JJti — 5) 2 + (t2 — 5) 2 < 4.5(ns), as shown in 
Fig. |j. Further, the Muon Counter hit informa- 
tion is used to identify di-muon events from other 
back-to-back two-prong final states, and this re- 



quires | cos < 0.65 due to the limited solid 
angle coverage. 

To separate electrons from muons and hadrons, 
the energies deposited by the two tracks in 
the Barrel Shower Counter (BSC) must sat- 
isfy ^ {E depl - If + {E dep2 - If < 0.65, where 

Edep = E b d " p is the normalized energy deposited. 
We also require |cos6* e | < 0.72, and because the 
Monte Carlo simulation does not model the en- 
ergy deposited well in the rib region of the BSC, 
an additional cut is applied on the z-coordinate 
of the first hit layer: 0.03 < \z sc \ < 0.85 or 
2 SC | > 0.95 m. 
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Figure 1. Times of flight for /i-pairs. The diag- 
onal bands are due to cosmic ray events. Good 
events, which cluster at the center of the plot are 
selected with the circle cut shown. 



A potential source of background is lepton- 
pairs coming from i/)(2S) — > XJ/ip, J/i/j — > . 
To eliminate this background, we make cuts on 
the track momenta, as shown in Fig. || p3| . 

Using the above criteria to select [i + and 
e + e~ events, we have compared various distri- 
butions with Monte Carlo distributions. Good 
agreement is found, as illustrated, for example, 
in Fig. |. 

For hadron event selection, we are guided by 
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our R-scan experience There is no par- 

ticular event topology to require; instead we 
make cuts to reject major backgrounds: cosmic 
rays, beam-associated background, two-photon 
processes (7*7*), mis-identified "hadron" events 
from QED processes of e + e~ — > I = e, 

H, t, and e + e~ — > 77 followed by 7 conver- 
sion, etc. Events with at least two well re- 
constructed charged tracks within |cos#| < 0.8 
are selected. The total energy deposited by an 
event in the BSC (E sum ) is required to be larger 
than 0.36-E{, eam , in order to suppress contamina- 
tion from two-photon processes and beam associ- 
ated background. Events with all tracks point- 
ing to the same hemisphere in the z direction 
are removed to suppress beam-associated back- 
ground. For two-prong events, two additional 
cuts are applied to eliminate possible lepton pair 
background. The number of photons must be 
greater than one, and the acollinearity between 
two charged tracks, otAcoU must be greater than 
10 degrees. A comparison between data and 
Monte Carlo samples satisfied the foregoing se- 
lection criteria is shown in Fig. |J. 

The background from r + r _ decay is difficult 
to distinguish from direct hadronic decay events, 
so the contribution from this source, N T + T -, is 
estimated using N TT — L ■ (e TT • &tt), where 
L is the integrated luminosity at each energy 
point, a TT the QED production cross section at 
this energy point, and e TT the acceptance of our 
hadron event selection criteria for t + t~ events. 
This is subtracted from the observed number 
of hadron events. A similar subtraction is per- 
formed for the other surviving backgrounds, such 
as e + e~ — > e + e~, 77, and the two-photon 

process (7*7*)- Therefore, the corrected number 
of hadron events, N£ bs , is 

s =N h - N TT - N ee - - A 77 - A 7 - 7 , 

where is the number of events that satisfy 
the hadron event selection cuts. The numbers of 
various backgrounds are given in Table |l|. 

For the selection of ip(2S) — > ir + ir~ J/ip events, 
we select a pair of low energy pions and deter- 
mine the mass recoiling against these two pions, 
m-recoii, which shows a strong J/ip peak, corre- 
sponding to the decay ip(2S) — > ir + ii~ J/ip. The 



Table 1 

Estimated numbers of hadron event back- 
grounds. 



E cm 
(GeV) 




Backg 


round events 




N TT 


Nee 


N„„ 


7V 77 


iV~»~. 
7 7 


3.6668 


20.79 


12.16 


0.22 


80 


1.45 


3.6719 


20.36 


11.91 


0.21 


78 


1.42 


3.6750 


20.15 


11.79 


0.21 


0.78 


1.41 


3.6781 


20.72 


12.12 


0.22 


0.80 


1.46 


3.6801 


20.58 


12.04 


0.21 


79 


1.45 


3.6809 


22.18 


12.98 


23 


85 


1.56 


3.6820 


20.80 


12.17 


0.22 


80 


1.46 


3.6828 


20.60 


12.05 


0.21 


79 


1.45 


3.6832 


20.66 


12.09 


0.22 


0.80 


1.46 


3.6836 


21.52 


12.59 


0.22 


0.83 


1.52 


3.6844 


20.61 


12.06 


0.22 


0.79 


1.45 


3.6850 


21.44 


12.54 


0.22 


0.83 


1.51 


3.6855 


19.98 


11.69 


0.21 


0.77 


1.41 


3.6863 


19.79 


11.58 


0.21 


0.76 


1.40 


3.6867 


19.92 


11.66 


0.21 


0.77 


1.41 


3.6875 


19.88 


11.63 


0.21 


0.77 


1.41 


3.6882 


20.42 


11.95 


0.21 


0.79 


1.44 


3.6886 


20.19 


11.81 


0.21 


0.78 


1.43 


3.6893 


25.36 


14.84 


0.26 


0.98 


1.79 


3.6908 


20.97 


12.27 


0.22 


0.81 


1.49 


3.6939 


23.11 


13.52 


0.24 


0.89 


1.64 


3.6979 


22.38 


13.09 


0.23 


0.86 


1.59 


3.7017 


21.89 


12.81 


0.23 


0.84 


1.56 


3.7068 


22.56 


13.20 


0.24 


0.87 


1.62 



mrecoii distribution is fitted with a signal shape 
plus polynomial background to obtain the num- 
ber of i>(2S) — > ir + ir~J/ip events at each en- 
ergy point. The very clean ip(2S) — > 7r + n~ J/ip, 
J/ip — » l + l~ channel is used to determine the sig- 
nal shape. The inclusive m r ecoii distribution for 
all 24 energy points combined is shown in Fig. |^. 
For more detail, see Rcf. [|l6). The numbers of 
events selected for the four final states at the 24 
energies are listed in Table |[ 

3. Acceptance 

The acceptance is the product of the trigger 
efficiency and the reconstruction-selection effi- 
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Table 2 

Numbers of events selected. 



Selected number of events 



(GeV) 


ATObS 

J 'ee 


pj-obs 
MM 


N obs 


ATObS 

7T+7T- J/lp 


3.6668 


1789 


110 


385 


7.2 


3.6719 


1752 


92 


389 


10.2 


3.6750 


1734 


124 


391 


22.3 


3.6781 


1783 


102 


437 


0.0 


3.6801 


1771 


66 


485 


12.8 


3.6809 


1909 


104 


588 


14.2 


3.6820 


1790 


91 


764 


75.0 


3.6828 


1773 


86 


1538 


181.0 


3.6832 


1778 


69 


2662 


426.1 


3.6836 


1852 


79 


4178 


668.0 


3.6844 


1774 


92 


8169 


1508.0 


3.6850 


1845 


128 


14684 


2660.8 


3.6855 


1720 


140 


15073 


2758.1 


3.6863 


1703 


191 


15621 


2902.6 


3.6867 


1715 


176 


14980 


2572.0 


3.6875 


1711 


161 


11658 


2201.8 


3.6882 


1758 


175 


7047 


1245.1 


3.6886 


1738 


125 


4960 


868.9 


3.6893 


2183 


188 


3964 


707.4 


3.6908 


1805 


113 


1897 


269.0 


3.6939 


1989 


97 


1318 


128.8 


3.6979 


1926 


132 


1029 


82.6 


3.7017 


1884 


118 


876 


77.2 


3.7068 


1942 


120 


729 


26.5 


sum 


43624 


2879 


113823 


19425.4 



ciency. The triggers are the same as those used 
in our R scan experiment [14 , The trigger ef- 
ficiencies, measured by comparing the responses 
to different trigger requirements in special runs 
taken at the J /ip resonance, are determined to 
be 1.000, 0.994 and 0.998 for e+e~,/i+/i- and 
hadronic events respectively, with an uncertainty 
of 0.005. 

Different generators are used to determine the 
reconstruction-selection efficiencies. For e + e~ 
and final states, the efficiencies for QED 

processes are determined in simulations with the 
BHABHA and MUPAIR generators @. The 
resonance e + e~ and fi + /i~ efficiencies are deter- 
mined using the generator V2LL, adapted from 



MUPAIR, with the initial state radiative correc- 
tions removed. For the hadronic processes, an 
event generator for charmonium inclusive decay 
flql is used to obtain the efficiency for the reso- 
nance portion, and the JETSET string fragmen- 
tation algorithm with parameters modified to fit 
the experimental data in the BEPC energy re- 
gion |L9| is used to compute the efficiency for 
the continuum portion. For the ir + Tr^J/'ip ac- 
ceptance, a phase space Monte Carlo program, 
modified to give the correct dipion mass and an- 
gular distributions pfj, is used. 



Table 3 

Acceptances for continuum, A c 
nance, A r . 



and the reso- 



final state 


e + e 




had 


7T+7T J/lj) 


A c (%) 


72.4 


37.1 


74.5 




SA C /A C (%) 


3.2 


5.0 


7.1 




A r {%) 


76.4 


41.9 


77.1 


43.4 


SA r /A r (%) 


8.5 


4.7 


2.2 


3.4 



The acceptances of the four final states for con- 
tinuum (A c ) and resonance (A r ) processes, to- 
gether with their relative errors, are listed in Ta- 
ble||. Here the acceptance for the e + e~ final state 
in the table applies to events within a restricted 
solid angle (| cos# e | < 0.72), whereas those of the 
hadron and dimuon final states cover all solid an- 
gles. The acceptance error includes the uncertain- 
ties estimated by varying selection cuts and us- 
ing different selection methods and Monte Carlo 
models. 

4. Fit of observed cross sections and re- 
sults 

The e + e~ — > hadrons, ir + ir~ J/ip, e + e~, and 
\x + pT events at the 24 scan points are fitted 
simultaneously to obtain the partial widths of 
ip(2S) to hadrons, n + ir~J/ip, and final 
states. The total width is assumed to be the sum 

+ r e + r r , 
r R = r„ = 



of four partial widths, T t = + 
and lepton universality is assumedy, 

*BES collaboration has measured the branching fraction 
of ip(2S) decay into t + t - . This value along with those 
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r r /0. 38847. A likelihood function is constructed 
22 1 . Two kinds of correlations are taken into con- 
sideration in the fitting formula: one is the corre- 
lation between different channels at the same en- 
ergy point, because the same luminosity is used to 
determine the cross sections at each energy; the 
other is the correlation between different energy 
points for each channel, due to the acceptance un- 
certainty being the same for all scan points. The 
M1NUIT package |^3| for maximization is used 
to give the best estimates for ip(2S) parameters 
and their uncertainties. The theoretical cross sec- 
tion used in the fit for the hadron channel uses 
a Breit-Wigner amplitude and a non-resonant 
direct-channel amplitude plus a J/ij) resonance 
"tail" cross section, as determined by a previous 
BES J/ip scan experiment [Q . The contributions 
from V'(25 l ) decay into e + e~, t + t~ final 

states mis-identified as hadron events have been 
taken into account by correcting the resonant 
hadron event acceptance^]. For the 7r + 7r~ J/tp fi- 
nal state, only a resonant Breit-Wigner amplitude 
is considered. In the fi + [i~ channel, the i/j(2S) 
resonant term, QED term, and their interference 
are included. The radiative corrections to these 
three processes are taken into account by the for- 
mulation of Refs. J2f| and [^6|. For the e + e~ final 
state, where the QED ^-channel photon exchange 
also contributes, a theoretical cross section in- 
cluding radiative corrections is derived using the 
method of Ref. [^7j . The effects on lepton-pair fi- 
nal state cross sections coming from vacuum po- 
larization are also taken into consideration p8[ ; 
while for hadronic final states, the vacuum polar- 
ization is absorbed into the definition of r ee , that 
is 

T exp (i;(2S) -> e+e-) = 



r (^(25) - e+e") 



ii-n(M| (2S) )r 



where T exp is the experimental width, Tq the low- 
est order in a width, and FJ the order a vacuum 



of the branching fractions in e+e — and /i + fj,~, satisfies 
the relation predicted by the sequential lepton hypothesis 
within errors |pH . 

^The contaminations from QED processes e + e _ — > l + l~ 
have been already subtracted, see section g. 



polarization |29[] . The theoretical cross sections 
are convoluted with the energy distribution of 
the colliding beams, which is treated as Gaus- 
sian. The following parameters are allowed to 
vary in the fit: the *p(2S) mass, M, the total 
width, r t , the partial widths, ^ w + w -j/^, and T/, 
the energy spread of the machine, and the non- 
resonant hadronic cross section. 

As the branching fraction of ip(2S) to e + e~ is 
small, the cross section for the e + e~ final state 
is dominated by the QED process. Therefore 
this channel is used to calculate the integrated 
luminosity at each energy point by an iterative 
method. First, all e + e~ events within |cos6*| < 
0.72 are taken as "Bhabha" events and used to 
calculate the integrated luminosity at each energy 
point. A maximum likelihood fit is performed to 
the observed cross sections for hadron, 7r + 7r~ J/tp, 
and muon pair final states, and a group of tp(2S) 
parameters is obtained with the assumption of 
e-fjL-r universality. Then, separating the e + e~ 
events into a QED part, a resonance part and 
their interference, the integrated luminosity for 
each energy point is recalculated using the QED 
part only, and the fitting procedures are redone 
to get new values for the ip(2S) parameters. The 
iterative process is repeated until the value of the 
integrated luminosity at each energy point is con- 
sistent for two successive iterations. 

The fitted curves are shown along with the scan 
points in Fig. ^. The fitted mass of the tp(2S) is 
corrected to the PDG value |50|. The errors in 
the other parameters caused by this correction are 
negligible. The fitted spread in the center-of-mass 
energy of the machine is (1.298 ± 0.007) MeV, 
in agreement with the expectation (~1.3 MeV). 
The resultant R ratio for the hadronic cross sec- 
tion near the ip(2S) resonance is 2.15 ± 0.16, 
which agrees well with the earlier BES R mea- 
surements p5| . 

The results of the fit for decay widths and 
branching fractions are given in Table ^, together 
with corresponding PDG ^] values for compar- 
ison. The errors are the sum in quadrature of 
statistical, fitting, and systematic uncertainties, 
including those from acceptance uncertainties, lu- 
minosity uncertainties, and a center-of-mass en- 
ergy uncertainty of 0.10 MeV [j32[. The system- 
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Table 4 



Results and comparison with the PDG2002 ]30| . 



Value 


BES 


PDG2002 


I t(kev) 1 


Z04 ± ZY 


o(JU ± zo 




(10.1%) 


/n o fry \ 

(8.3%) 


i h (kev j 


oc;o _i_ OA 
ZOO it ZD 

(10.1%) 




r 7r + T -j/v,(keV) t 


85.4 ±8.7 
(10.1%) 




r ; (keV)f 


2.44 ±0.21 


2.19 ±0.15 




(8.8%) 


(6.8%)* 




97.79 ±0.15 


98.10 ±0.30 




(0.16%) 


(0.31 %) 




32.3 ± 1.4 
(4.4%) 


30.5 ± 1.6 
(5.2 %) 




0.93 ±0.08 


0.73 ±0.04 




(8.5%) 


(5.5%)* 



Note : f indicates directly fitted values while others are 
derived quantities; J indicates the PDG values for T e 
and B e - The numbers in parenthesis denote the relative 
errors. 



atic error for T t , T^+^-j/^ and Ti is (2.24 — 4.7) 
%, 3.2 % and (6.4 - 7.0) %, respectively. The 
correlation coefficients between directly fitted pa- 
rameters obtained from the fitting are utilized in 
the calculation of the errors of indirectly deter- 
mined parameters such as I\ and branching frac- 
tions. For Bh, B^+k- jMj and Bi, the uncertainty 
related to the common error from the luminosity 
measurement cancels out. 

As a check of the fitting procedure, B^+^-j/.^ 
was also determined by a simpler approach. In 
this approach, the distribution of N% bs versus 
energy was fit with the shape determined from 
^V°+ S 7r - versus scan energy plus a polynomial 
to represent the continuum process to determine 
the number of hadrons coming from ip(2S) de- 
cays. Using the ratio of the total number of 
iV°+ 7r _ events corrected by their detection effi- 
ciency and the number of 4>{2S) hadronic decays 
corrected by their detection efficiency, we directly 
determine B n + n -j/^,, which agrees very well with 
the result from the full fitting procedure. 

The assumption that lepton pairs couple to the 
ip(2S) only via an intermediate photon [ |3l| im- 



plies the existence of the decay ip(2S) — > 7 — > 
hadrons with a branching fraction: 

T y h/Tt = RTjTt = 0.0199 ± 0.0019, 

which corresponds to a width Y 1 h of 5.26 ± 0.32 
keV. 

The ip(2S) lepton width T ; = (2.44 ± 0.21) 
keV obtained in this measurement agrees with 
the BES previous value of T e = (2.07 ± 0.32) 
keV J8|, within the error. In addition, this is 
the first direct measurement to the decay width 
of ip(2S) — » 7r + 7r~ J ftp, and the precision of 
Bn+n- j/ii> is much better than previous measure- 
ments and the current PDG value. 
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Figure 3. Normalized momentum distributions 
for (a) e + e~ events and (b) fi + (J.~ events. (His- 
togram for Monte Carlo events and dots with er- 
ror bar for data) 



Figure 2. Momentum distributions for (a) e + e~ 
events and (b) events. The solid line indi- 

cates the cuts applied to remove background from 



9 




o.o-i 
a. acts 



:== O.D3 

a 

£> O.O&S 

o.o-is 
a. a i 





w 




(b) 
















_ \ 

























-It) SO SO lOO I^O I -It) IOO ISO 



it ) 



Figure 4. Comparison between Monte Carlo 
events (histogram) and data (dots with error bar) 
for hadron events, (a) Esum (b) ctAcoi (c) cos 9 
of charged track. (The arrow in (a) indicates the 
cut applied.) 
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Figure 5. Mass recoiling against 7r + 7r~ for all 24 
energy points combined. The J/ip peak is due to 
the signal tp(2S) — > -k + 7t~ J/ip, J/ip -^anything 
events. 
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Figure 6. The cross section for (a) e+e~ — > 
hadrons, (b) e + e~ — > ir + ir~ J/ip, and (c) e + e~ — > 
versus center-of-mass energy. The solid 
curves represent the results of the fit to the data. 



